Force-free state in a superconducting single crystal and angle-dependent vortex helical instability

Superconducting 2H−NbSe2 single crystals show intrinsic low pinning values. Therefore, they are ideal materials with which to explore fundamental properties of vortices. (V, I) characteristics are the experimental data we have used to investigate the dissipation mechanisms in a rectangular-shaped 2H−NbSe2 single crystal. Particularly, we have studied dissipation behavior with magnetic fields applied in the plane of the crystal and parallel to the injected currents, i.e., in the force-free state where the vortex helical instability governs the vortex dynamics. In this regime, the data follow the elliptic critical state model and the voltage dissipation shows an exponential dependence, V∝eα(I−IC∥), IC∥ being the critical current in the force-free configuration and α a linear temperature-dependent parameter. Moreover, this exponential dependence can be observed for in-plane applied magnetic fields up to 40° off the current direction, which implies that the vortex helical instability plays a role in dissipation even out of the force-free configuration

Ultra-broad spectral photo-response in FePS3 air-stable devices

Van der Waals materials with narrow energy gaps and efficient response over a broadband optical spectral range are key to widen the energy window of nanoscale optoelectronic devices. Here, we characterize FePS3 as an appealing narrow-gap p-type semiconductor with an efficient broadband photo-response, a high refractive index, and a remarkable resilience against air and light exposure. To enable fast prototyping, we provide a straightforward guideline to determine the thickness of few-layered FePS3 nanosheets extracted from the optical transmission characteristics of several flakes. The analysis of the electrical photo-response of FePS3 devices as a function of the excitation energy confirms a narrow gap suitable for near IR detection (1.23 eV) and, more interestingly, reveals a broad spectral responsivity up to the ultraviolet region. The experimental estimate for the gap energy is corroborated by ab-initio calculations. An analysis of photocurrent as a function of gate voltage and incident power reveals a photo-response dominated by photogating effects. Finally, aging studies of FePS3 nanosheets under ambient conditions show a limited reactivity of the outermost layers of flakes in long exposures to air

Computational and in vitro Pharmacodynamics Characterization of 1A-116 Rac1 Inhibitor: Relevance of Trp56 in Its Biological Activity

In the last years, the development of new drugs in oncology has evolved notably. In particular, drug development has shifted from empirical screening of active cytotoxic compounds to molecularly targeted drugs blocking specific biologic pathways that drive cancer progression and metastasis. Using a rational design approach, our group has developed 1A-116 as a promising Rac1 inhibitor, with antitumoral and antimetastatic effects in several types of cancer. Rac1 is over activated in a wide range of tumor types and and it is one of the most studied proteins of the Rho GTPase family. Its role in actin cytoskeleton reorganization has effects on endocytosis, vesicular trafficking, cell cycle progression and cellular migration. In this context, the regulatory activity of Rac1 affects several key processes in the course of the cancer including invasion and metastasis. The purpose of this preclinical study was to focus on the mode of action of 1A-116, conducting an interdisciplinary approach with in silico bioinformatics tools and in vitro assays. Here, we demonstrate that the tryptophan 56 residue is necessary for the inhibitory effects of 1A-116 since this compound interferes with protein-protein interactions (PPI) of Rac1GTPase involving several GEF activators. 1A116 is also able to inhibit the oncogenic Rac1P29S mutant protein, one of the oncogenic drivers found in sun-exposed melanoma. It also inhibits numerous Rac1-regulated cellular processes such as membrane ruffling and lamellipodia formation. These results deepen our knowledge of 1A-116 inhibition of Rac1 and its biological impact on cancer progression. They also represent a good example of how in silico analyses represent a valuable approach for drug development.Fil: González, Nazareno. Universidad Nacional de Quilmes. Departamento de Ciencia y Tecnología. Laboratorio de Oncología Molecular; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; ArgentinaFil: Cardama, Georgina Alexandra. Universidad Nacional de Quilmes. Departamento de Ciencia y Tecnología. Laboratorio de Oncología Molecular; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; ArgentinaFil: Chinestrad, Patricio Manuel. Universidad Nacional de Quilmes. Departamento de Ciencia y Tecnología. Laboratorio de Farmacología Molecular; ArgentinaFil: Robles Valero, Javier. Universidad de Salamanca; España. Consejo Superior de Investigaciones Científicas; EspañaFil: Rodríguez Fdez, Sonia. Universidad de Salamanca; España. Consejo Superior de Investigaciones Científicas; EspañaFil: Lorenzo Martín, L. Francisco. Universidad de Salamanca; España. Consejo Superior de Investigaciones Científicas; EspañaFil: Bustelo, Xosé R.. Universidad de Salamanca; España. Consejo Superior de Investigaciones Científicas; EspañaFil: Lorenzano Menna, Pablo. Universidad Nacional de Quilmes. Departamento de Ciencia y Tecnología. Laboratorio de Farmacología Molecular; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; ArgentinaFil: Gomez, Daniel Eduardo. Universidad Nacional de Quilmes. Departamento de Ciencia y Tecnología. Laboratorio de Oncología Molecular; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; Argentin

Analysis of the Legionella longbeachae Genome and Transcriptome Uncovers Unique Strategies to Cause Legionnaires' Disease

Legionella pneumophila and L. longbeachae are two species of a large genus of bacteria that are ubiquitous in nature. L. pneumophila is mainly found in natural and artificial water circuits while L. longbeachae is mainly present in soil. Under the appropriate conditions both species are human pathogens, capable of causing a severe form of pneumonia termed Legionnaires' disease. Here we report the sequencing and analysis of four L. longbeachae genomes, one complete genome sequence of L. longbeachae strain NSW150 serogroup (Sg) 1, and three draft genome sequences another belonging to Sg1 and two to Sg2. The genome organization and gene content of the four L. longbeachae genomes are highly conserved, indicating strong pressure for niche adaptation. Analysis and comparison of L. longbeachae strain NSW150 with L. pneumophila revealed common but also unexpected features specific to this pathogen. The interaction with host cells shows distinct features from L. pneumophila, as L. longbeachae possesses a unique repertoire of putative Dot/Icm type IV secretion system substrates, eukaryotic-like and eukaryotic domain proteins, and encodes additional secretion systems. However, analysis of the ability of a dotA mutant of L. longbeachae NSW150 to replicate in the Acanthamoeba castellanii and in a mouse lung infection model showed that the Dot/Icm type IV secretion system is also essential for the virulence of L. longbeachae. In contrast to L. pneumophila, L. longbeachae does not encode flagella, thereby providing a possible explanation for differences in mouse susceptibility to infection between the two pathogens. Furthermore, transcriptome analysis revealed that L. longbeachae has a less pronounced biphasic life cycle as compared to L. pneumophila, and genome analysis and electron microscopy suggested that L. longbeachae is encapsulated. These species-specific differences may account for the different environmental niches and disease epidemiology of these two Legionella species

More than 18,000 effectors in the Legionella genus genome provide multiple, independent combinations for replication in human cells.

The genus Legionella comprises 65 species, among which Legionella pneumophila is a human pathogen causing severe pneumonia. To understand the evolution of an environmental to an accidental human pathogen, we have functionally analyzed 80 Legionella genomes spanning 58 species. Uniquely, an immense repository of 18,000 secreted proteins encoding 137 different eukaryotic-like domains and over 200 eukaryotic-like proteins is paired with a highly conserved type IV secretion system (T4SS). Specifically, we show that eukaryotic Rho- and Rab-GTPase domains are found nearly exclusively in eukaryotes and Legionella Translocation assays for selected Rab-GTPase proteins revealed that they are indeed T4SS secreted substrates. Furthermore, F-box, U-box, and SET domains were present in >70% of all species, suggesting that manipulation of host signal transduction, protein turnover, and chromatin modification pathways are fundamental intracellular replication strategies for legionellae. In contrast, the Sec-7 domain was restricted to L. pneumophila and seven other species, indicating effector repertoire tailoring within different amoebae. Functional screening of 47 species revealed 60% were competent for intracellular replication in THP-1 cells, but interestingly, this phenotype was associated with diverse effector assemblages. These data, combined with evolutionary analysis, indicate that the capacity to infect eukaryotic cells has been acquired independently many times within the genus and that a highly conserved yet versatile T4SS secretes an exceptional number of different proteins shaped by interdomain gene transfer. Furthermore, we revealed the surprising extent to which legionellae have coopted genes and thus cellular functions from their eukaryotic hosts, providing an understanding of how dynamic reshuffling and gene acquisition have led to the emergence of major human pathogens

Accurate reconstruction of insertion-deletion histories by statistical phylogenetics

The Multiple Sequence Alignment (MSA) is a computational abstraction that represents a partial summary either of indel history, or of structural similarity. Taking the former view (indel history), it is possible to use formal automata theory to generalize the phylogenetic likelihood framework for finite substitution models (Dayhoff's probability matrices and Felsenstein's pruning algorithm) to arbitrary-length sequences. In this paper, we report results of a simulation-based benchmark of several methods for reconstruction of indel history. The methods tested include a relatively new algorithm for statistical marginalization of MSAs that sums over a stochastically-sampled ensemble of the most probable evolutionary histories. For mammalian evolutionary parameters on several different trees, the single most likely history sampled by our algorithm appears less biased than histories reconstructed by other MSA methods. The algorithm can also be used for alignment-free inference, where the MSA is explicitly summed out of the analysis. As an illustration of our method, we discuss reconstruction of the evolutionary histories of human protein-coding genes.Comment: 28 pages, 15 figures. arXiv admin note: text overlap with arXiv:1103.434

Activation of Type 1 Cannabinoid Receptor (CB1R) promotes neurogenesis in murine subventricular zone cell cultures

The endocannabinoid system has been implicated in the modulation of adult neurogenesis. Here, we describe the effect of type 1 cannabinoid receptor (CB1R) activation on self-renewal, proliferation and neuronal differentiation in mouse neonatal subventricular zone (SVZ) stem/progenitor cell cultures. Expression of CB1R was detected in SVZ-derived immature cells (Nestin-positive), neurons and astrocytes. Stimulation of the CB1R by (R)-(+)-Methanandamide (R-m-AEA) increased self-renewal of SVZ cells, as assessed by counting the number of secondary neurospheres and the number of Sox2+/+ cell pairs, an effect blocked by Notch pathway inhibition. Moreover, R-m-AEA treatment for 48 h, increased proliferation as assessed by BrdU incorporation assay, an effect mediated by activation of MAPK-ERK and AKT pathways. Surprisingly, stimulation of CB1R by R-m-AEA also promoted neuronal differentiation (without affecting glial differentiation), at 7 days, as shown by counting the number of NeuN-positive neurons in the cultures. Moreover, by monitoring intracellular calcium concentrations ([Ca2+](i)) in single cells following KCl and histamine stimuli, a method that allows the functional evaluation of neuronal differentiation, we observed an increase in neuronal-like cells. This proneurogenic effect was blocked when SVZ cells were co-incubated with R-m-AEA and the CB1R antagonist AM 251, for 7 days, thus indicating that this effect involves CB1R activation. In accordance with an effect on neuronal differentiation and maturation, R-m-AEA also increased neurite growth, as evaluated by quantifying and measuring the number of MAP2-positive processes. Taken together, these results demonstrate that CB1R activation induces proliferation, self-renewal and neuronal differentiation from mouse neonatal SVZ cell cultures.Fundacao para a Ciencia e a Tecnologia - Portugal [POCTI/SAU-NEU/68465/2006, PTDC/SAU-NEU/104415/2008, PTDC/SAU-NEU/101783/2008, POCTI/SAU-NEU/110838/2009]; Fundacao Calouste Gulbenkian [96542]; Fundacao para a Ciencia e Tecnologiainfo:eu-repo/semantics/publishedVersio

Metabolic Networks of Sodalis glossinidius: A Systems Biology Approach to Reductive Evolution

Background: Genome reduction is a common evolutionary process affecting bacterial lineages that establish symbiotic or pathogenic associations with eukaryotic hosts. Such associations yield highly reduced genomes with greatly streamlined metabolic abilities shaped by the type of ecological association with the host. Sodalis glossinidius, the secondary endosymbiont of tsetse flies, represents one of the few complete genomes available of a bacterium at the initial stages of this process. In the present study, genome reduction is studied from a systems biology perspective through the reconstruction and functional analysis of genome-scale metabolic networks of S. glossinidius. Results: The functional profile of ancestral and extant metabolic networks sheds light on the evolutionary events underlying transition to a host-dependent lifestyle. Meanwhile, reductive evolution simulations on the extant metabolic network can predict possible future evolution of S. glossinidius in the context of genome reduction. Finally, knockout simulations in different metabolic systems reveal a gradual decrease in network robustness to different mutational events for bacterial endosymbionts at different stages of the symbiotic association. Conclusions: Stoichiometric analysis reveals few gene inactivation events whose effects on the functionality of S. glossinidius metabolic systems are drastic enough to account for the ecological transition from a free-living to hostdependent lifestyle. The decrease in network robustness across different metabolic systems may be associated with th

Chance and necessity in the genome evolution of endosymbiotic bacteria of insects

[EN] An open question in evolutionary biology is how does the selection¿drift balance determine the fates of biological interactions. We searched for signatures of selection and drift in genomes of five endosymbiotic bacterial groups known to evolve under strong genetic drift. Although most genes in endosymbiotic bacteria showed evidence of relaxed purifying selection, many genes in these bacteria exhibited stronger selective constraints than their orthologs in free-living bacterial relatives. Remarkably, most of these highly constrained genes had no role in the host¿symbiont interactions but were involved in either buffering the deleterious consequences of drift or other host-unrelated functions, suggesting that they have either acquired new roles or their role became more central in endosymbiotic bacteria. Experimental evolution of Escherichia coli under strong genetic drift revealed remarkable similarities in the mutational spectrum, genome reduction patterns and gene losses to endosymbiotic bacteria of insects. Interestingly, the transcriptome of the experimentally evolved lines showed a generalized deregulation of the genome that affected genes encoding proteins involved in mutational buffering, regulation and amino acid biosynthesis, patterns identical to those found in endosymbiotic bacteria. Our results indicate that drift has shaped endosymbiotic associations through a change in the functional landscape of bacterial genes and that the host had only a small role in such a shiftThis work was supported by Science Foundation Ireland (12/IP/1637) and grants from the Spanish Ministerio de Economia y Competitividad (MINECO-FEDER; BFU2012-36346 and BFU2015-66073-P) to MAF. DAP and CT were supported by Juan de la Cierva fellowships from MINECO (references: JCI-2011-11089 and JCA-2012-14056, respectively). DAP is supported by funds from the University of Nevada, Reno, NV, USA.Sabater-Muñoz, B.; Toft, C.; Alvarez-Ponce, D.; Fares Riaño, MA. (2017). Chance and necessity in the genome evolution of endosymbiotic bacteria of insects. The ISME Journal. 11(6):1291-1304. https://doi.org/10.1038/ismej.2017.18S12911304116Aguilar-Rodriguez J, Sabater-Munoz B, Montagud-Martinez R, Berlanga V, Alvarez-Ponce D, Wagner A et al. (2016). The molecular chaperone DnaK is a source of mutational robustness. Genome Biol Evol 8: 2979–2991.Alvarez-Ponce D, Sabater-Munoz B, Toft C, Ruiz-Gonzalez MX, Fares MA . (2016). Essentiality is a strong determinant of protein rates of evolution during mutation accumulation experiments in Escherichia coli. Genome Biol Evol 8: 2914–2927.Anders S, Huber W . (2010). Differential expression analysis for sequence count data. Genome Biol 11: R106.Archibald J . (2014) One Plus One Equals One: Symbiosis and the Evolution of Complex Life. Oxford University Press: Oxford, UK.Aussel L, Loiseau L, Hajj Chehade M, Pocachard B, Fontecave M, Pierrel F et al. (2014). ubiJ, a new gene required for aerobic growth and proliferation in macrophage, is involved in coenzyme Q biosynthesis in Escherichia coli and Salmonella enterica serovar Typhimurium. J Bacteriol 196: 70–79.Baumann P, Baumann L, Clark MA . (1996). Levels of Buchnera aphidicola chaperonin groEL during growth of the aphid Schizaphis graminum. Curr Microbiol 32: 7.Benjamini Y, Yekutieli Y . (2005). False discovery rate controlling confidence intervals for selected parameters. J Am Stat Assoc 100: 10.Bennett GM, Moran NA . (2015). Heritable symbiosis: the advantages and perils of an evolutionary rabbit hole. Proc Natl Acad Sci USA 112: 10169–10176.Bermingham J, Rabatel A, Calevro F, Vinuelas J, Febvay G, Charles H et al. (2009). Impact of host developmental age on the transcriptome of the symbiotic bacterium Buchnera aphidicola in the pea aphid (Acyrthosiphon pisum. Appl Environ Microbiol 75: 7294–7297.Bogumil D, Dagan T . (2010). Chaperonin-dependent accelerated substitution rates in prokaryotes. Genome Biol Evol 2: 602–608.Carbon S, Ireland A, Mungall CJ, Shu S, Marshall B, Lewis S et al. (2009). AmiGO: online access to ontology and annotation data. Bioinformatics 25: 288–289.Chen Z, Wang Y, Li Y, Li Y, Fu N, Ye J et al. (2012). Esre: a novel essential non-coding RNA in Escherichia coli. FEBS Lett 586: 1195–1200.Clark JW, Hossain S, Burnside CA, Kambhampati S . (2001). Coevolution between a cockroach and its bacterial endosymbiont: a biogeographical perspective. Proc Biol Sci 268: 393–398.Dale C, Wang B, Moran N, Ochman H . (2003). Loss of DNA recombinational repair enzymes in the initial stages of genome degeneration. Mol Biol Evol 20: 1188–1194.Deatherage DE, Barrick JE . (2014). Identification of mutations in laboratory-evolved microbes from next-generation sequencing data using breseq. Methods Mol Biol 1151: 165–188.Douglas AE . (2003). The nutritional physiology of aphids. Adv Insect Physiol 31: 68.Fares MA, Barrio E, Sabater-Munoz B, Moya A . (2002a). The evolution of the heat-shock protein GroEL from Buchnera, the primary endosymbiont of aphids, is governed by positive selection. Mol Biol Evol 19: 1162–1170.Fares MA, Ruiz-Gonzalez MX, Moya A, Elena SF, Barrio E . (2002b). Endosymbiotic bacteria: groEL buffers against deleterious mutations. Nature 417: 398.Gancedo C, Flores CL, Gancedo JM . (2016). The expanding landscape of moonlighting proteins in yeasts. Microbiol Mol Biol Rev 80: 765–777.Gerardo NM, Altincicek B, Anselme C, Atamian H, Barribeau SM, de Vos M et al. (2010). Immunity and other defenses in pea aphids, Acyrthosiphon pisum. Genome Biol 11: R21.Gomez-Valero L, Latorre A, Silva FJ . (2004). The evolutionary fate of nonfunctional DNA in the bacterial endosymbiont Buchnera aphidicola. Mol Biol Evol 21: 2172–2181.Gomez-Valero L, Silva FJ, Christophe Simon J, Latorre A . (2007). Genome reduction of the aphid endosymbiont Buchnera aphidicola in a recent evolutionary time scale. Gene 389: 87–95.Gonzalez-Domenech CM, Belda E, Patino-Navarrete R, Moya A, Pereto J, Latorre A . (2012). Metabolic stasis in an ancient symbiosis: genome-scale metabolic networks from two Blattabacterium cuenoti strains, primary endosymbionts of cockroaches. BMC Microbiol 12 (Suppl 1): S5.Hansen AK, Moran NA . (2011). Aphid genome expression reveals host-symbiont cooperation in the production of amino acids. Proc Natl Acad Sci USA 108: 2849–2854.Hansen AK, Moran NA . (2014). The impact of microbial symbionts on host plant utilization by herbivorous insects. Mol Ecol 23: 1473–1496.Henderson B, Fares MA, Lund PA . (2013). Chaperonin 60: a paradoxical, evolutionarily conserved protein family with multiple moonlighting functions. Biol Rev Camb Philos Soc 88: 955–987.Humphreys NJ, Douglas AE . (1997). Partitioning of symbiotic bacteria between generations of an insect: a quantitative study of a Buchnera sp. in the pea aphid (Acyrthosiphon pisum reared at different temperatures. Appl Environ Microbiol 63: 3294–3296.International Aphid Genomics Consortium. (2010). Genome sequence of the pea aphid Acyrthosiphon pisum. PLoS Biol 8: e1000313.Kadibalban AS, Bogumil D, Landan G, Dagan T . (2016). DnaK-dependent accelerated evolutionary rate in prokaryotes. Genome Biol Evol 8: 1590–1599.Katoh K, Standley DM . (2013). MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol 30: 772–780.Kelkar YD, Ochman H . (2013). Genome reduction promotes increase in protein functional complexity in bacteria. Genetics 193: 303–307.Koga R, Meng XY, Tsuchida T, Fukatsu T . (2012). Cellular mechanism for selective vertical transmission of an obligate insect symbiont at the bacteriocyte-embryo interface. Proc Natl Acad Sci USA 109: E1230–E1237.Kuo CH, Moran NA, Ochman H . (2009). The consequences of genetic drift for bacterial genome complexity. Genome Res 19: 1450–1454.Kuo CH, Ochman H . (2009). Deletional bias across the three domains of life. Genome Biol Evol 1: 145–152.Law R, Lewis DH . (1983). Biotic environments and the maintenance of sex-some evidence from mutualistic symbioses. Biol J Linnean Soc 20: 28.Liu XD, Xie L, Wei Y, Zhou X, Jia B, Liu J et al. (2014). Abiotic stress resistance, a novel moonlighting function of ribosomal protein RPL44 in the halophilic fungus Aspergillus glaucus. Appl Environ Microbiol 80: 4294–4300.Lohse M, Bolger AM, Nagel A, Fernie AR, Lunn JE, Stitt M et al. (2012). RobiNA: a user-friendly, integrated software solution for RNA-Seq-based transcriptomics. Nucleic Acids Res 40: W622–W627.Macdonald SJ, Lin GG, Russell CW, Thomas GH, Douglas AE . (2012). The central role of the host cell in symbiotic nitrogen metabolism. Proc Biol Sci 279: 2965–2973.McClure R, Balasubramanian D, Sun Y, Bobrovskyy M, Sumby P, Genco CA et al. (2013). Computational analysis of bacterial RNA-Seq data. Nucleic Acids Res 41: e140.McCutcheon JP, Moran NA . (2012). Extreme genome reduction in symbiotic bacteria. Nat Rev Microbiol 10: 13–26.McFall-Ngai M, Hadfield MG, Bosch TC, Carey HV, Domazet-Loso T, Douglas AE et al. (2013). Animals in a bacterial world, a new imperative for the life sciences. Proc Natl Acad Sci USA 110: 3229–3236.Mira A, Ochman H, Moran NA . (2001). Deletional bias and the evolution of bacterial genomes. Trends Genet 17: 589–596.Moran NA . (1996). Accelerated evolution and Muller's rachet in endosymbiotic bacteria. Proc Natl Acad Sci USA 93: 2873–2878.Moran NA, Dunbar HE, Wilcox JL . (2005). Regulation of transcription in a reduced bacterial genome: nutrient-provisioning genes of the obligate symbiont Buchnera aphidicola. J Bacteriol 187: 4229–4237.Moran NA, McCutcheon JP, Nakabachi A . (2008). Genomics and evolution of heritable bacterial symbionts. Annu Rev Genet 42: 165–190.Moran NA, McLaughlin HJ, Sorek R . (2009). The dynamics and time scale of ongoing genomic erosion in symbiotic bacteria. Science 323: 379–382.Nakabachi A, Ishida K, Hongoh Y, Ohkuma M, Miyagishima SY . (2014). Aphid gene of bacterial origin encodes a protein transported to an obligate endosymbiont. Curr Biol 24: R640–R641.Nilsson AI, Koskiniemi S, Eriksson S, Kugelberg E, Hinton JC, Andersson DI . (2005). Bacterial genome size reduction by experimental evolution. Proc Natl Acad Sci USA 102: 12112–12116.Patino-Navarrete R, Moya A, Latorre A, Pereto J . (2013). Comparative genomics of Blattabacterium cuenoti: the frozen legacy of an ancient endosymbiont genome. Genome Biol Evol 5: 351–361.Pettersson ME, Berg OG . (2007). Muller's ratchet in symbiont populations. Genetica 130: 199–211.Price DR, Feng H, Baker JD, Bavan S, Luetje CW, Wilson AC . (2014). Aphid amino acid transporter regulates glutamine supply to intracellular bacterial symbionts. Proc Natl Acad Sci USA 111: 320–325.Reyes-Prieto M, Vargas-Chavez C, Latorre A, Moya A . (2015). SymbioGenomesDB: a database for the integration and access to knowledge on host-symbiont relationships. Database 2015: bav109 (1–8).Robinson MD, McCarthy DJ, Smyth GK . (2010). edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26: 139–140.Sabater-Muñoz B, Prats-Escriche M, Montagud-Martinez R, Lopez-Cerdan A, Toft C, Aguilar-Rodriguez J et al. (2015). Fitness trade-offs determine the role of the molecular chaperonin groel in buffering mutations. Mol Biol Evol 32: 2681–2693.Schlicker A, Domingues FS, Rahnenfuhrer J, Lengauer T . (2006). A new measure for functional similarity of gene products based on Gene Ontology. BMC Bioinformatics 7: 302.Shigenobu S, Watanabe H, Hattori M, Sakaki Y, Ishikawa H . (2000). Genome sequence of the endocellular bacterial symbiont of aphids Buchnera sp. APS. Nature 407: 81–86.Supek F, Bosnjak M, Skunca N, Smuc T . (2011). REVIGO summarizes and visualizes long lists of gene ontology terms. PLoS ONE 6: e21800.Tamas I, Klasson L, Canback B, Naslund AK, Eriksson AS, Wernegreen JJ et al. (2002). 50 million years of genomic stasis in endosymbiotic bacteria. Science 296: 2376–2379.Toft C, Fares MA . (2008). The evolution of the flagellar assembly pathway in endosymbiotic bacterial genomes. Mol Biol Evol 25: 2069–2076.van Ham RC, Kamerbeek J, Palacios C, Rausell C, Abascal F, Bastolla U et al. (2003). Reductive genome evolution in Buchnera aphidicola. Proc Natl Acad Sci USA 100: 581–586.Wernegreen JJ . (2002). Genome evolution in bacterial endosymbionts of insects. Nat Rev Genet 3: 850–861.Wernegreen JJ . (2011). Reduced selective constraint in endosymbionts: elevation in radical amino acid replacements occurs genome-wide. PLoS One 6: e28905.Williams TA, Fares MA . (2010). The effect of chaperonin buffering on protein evolution. Genome Biol Evol 2: 609–619.Yang Z . (2007). PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol 24: 1586–1591